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We have precisely determined the ground state phase diagram of the quantum spin-1 bihnear- 
biquadratic Heisenberg model on the honeycomb lattice using the tensor renormalization group 
method. We find that the ferromagnetic, ferroquadrupolar, and a large part of the antiferromagnetic 
phases are stable against quantum fluctuations. However, around the phase where the ground state 
is antiferro-quadrupolar ordered in the classical limit, quantum fluctuations suppress completely all 
magnetic orders, leading to a plaquette order phase which breaks the lattice symmetry but preserves 
the spin SU(2) symmetry. On the evidence of our numerical results, the quantum phase transition 
between the antiferromagnetic phase and the plaquette phase is found to be either a direct second 
order or a very weak first order transition. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg 



I. INTRODUCTION 

In a quantum spin system, the spin order of the classi- 
cal ground state can be melted by a quantum fluctuation 
at zero temperature^, and the resulting so-called quan- 
tum spin liquid has been suggested as a possible parent 
state of high temperature superconductivity upon elec- 
tron or hole doping^ . The quantum fluctuation is usually 
enhanced with small spin, low dimensionality, and geo- 
metric frustration. Indeed, so far almost all the candi- 
dates of quantum spin liquid discovered or proposed are 
effective spin- 1/2 systems with or without charge fluc- 
tuations on triangular lattice^—, on Kagome lattice^^^, 
or on honeycomb lattice^^-— . Besides the exotic quan- 
tum spin liquid phases, strong quantum fluctuations in 
spin- 1/2 systems can also lead to highly unconventional 
quantum critical points. For example, it has been pro- 
posed that a generic direct second order quantum phase 
transition between a magnetic ordered phase and a para- 
magnetic phase with broken lattice symmetry can exist 
in spin- 1/2 quantum magnetsi^iii^. Such a transition is 
forbidden by the classic Landau-Ginzburg- Wilson-Fisher 
paradigm, and is called the deconfined quantum critical 
point, as it is described by fractionalized quantities in- 
stead of physical order parameters. In the last few years 
this theoretical proposal has gained strong numerical ev- 
idence by quantum Monte Carlo simulation on spin- 1/2 
models with both nearest neighbor Heisenberg coupling 
and four-spin interactions^Sri^^. 

In this paper we will address the question: Can quan- 
tum fluctuation lead to such exotic phases and phase 
transitions in systems with larger spins in two dimen- 
sions without geometric frustration? Theoretically, this 
question is highly nontrivial as the Affleck-Kennedy-Lieb- 
Tasaki type of valence bond state (VBS) only exists for 
spin-1 systems in one dimension. Recent experiments on 



the spin-1 magnet Ba3NiSb209 suggested that a highly 
nontrivial quantum disordered ground state of a two di- 
mensional spin-1 system is indeed possible^^. 

By using the state of the art tensor renormalization 
group method^^^^, we present strong numerical evi- 
dences for the quantum fluctuation driven exotic physics 
in the spin-1 bilinear-biquadratic Heisenberg model on 
the honeycomb lattice. An interesting but surprising re- 
sult we find is that all dipole and quadruple magnetic 
orders vanish in a phase where the ground state is stag- 
gered quadrupolar ordered in the classical limit, instead 
the system develops a translation symmetry breaking pla- 
quette order. Moreover, we demonstrate that the tran- 
sition between the plaquette and antiferromagnetic AF 
order is either a direct second order transition or a very 
weak first order transition. If it is indeed a direct sec- 
ond order transition, then it is most likely a deconfined 
quantum critical point, which is the first example of de- 
confined quantum critical point in spin-1 system. 

The spin-1 bilinear-biquadratic Heisenberg model 
reads 



H 



S, + (sin^)( -S,)^ 



(1) 



The honeycomb lattice has the smallest coordination 
number in two dimensions and the effect of quantum fluc- 
tuation is the strongest. This model contains a number 
of special points. The point ^ = is the conventional 
SU(2) AF Heisenberg model. When (9 = 7r/4, 57r/4, 
or ±7r/2, the Hamiltonian is SU(3) invariant, possess- 
ing a symmetry higher than the spin SU(2) symmetry. 
At 0± = ±arctan2, the Hamiltonian can be expressed 
purely using the quadrupolar tensor operator 
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Like the ferromagnetic (FM) spin operator, the uniform 
quadrupolar operator, Q = Qi, commutes with this 
Hamiltonian. However, the staggered quadruple opera- 
tor, Qg = Xli(~)*Q*5 does not commute with the Hamil- 
tonian. 

The two terms in Eq. ^ introduce competition be- 
tween different kinds of magnetic orders. The first term 
favors the conventional ferromagnetic or antiferromag- 
netic order, while the second term favors a ferro- or 
antiferro-quadrupolar order. This competition causes 
a strong quantum fluctuation, especially in the regime 
sin^ > where the Marshall sign rule is not applica- 
ble to the ground state wave function and the quantum 
Monte Carlo suffers the minus-sign problem. 

Aspects of the spin-1 bilinear biquadratic model have 
been explored previously in the literature. In one dimen- 
sion, the ground state phase diagram has been character- 
ized by numerical density matrix renormalization group 
method. For — 7r/4 < < 7r/4, the model gives rise to the 
Haldane spin gapped phase, while the ground state for 
7r/4 < 6 < 7r/2 corresponds to a quantum critical phase 
with power-law spin and quadrupolar correlations^^'^^. 

In two dimensions, the ground state phase diagram has 
not been firmly established. In the classical limit, this 
model possesses four phases^^'^^, as depicted by the inner 
circle of Fig. [TJ In the lower half plane of ^, the quan- 
tum Monte Carlo simulation^^ and other calculations^^ 
confirmed the classical phase diagram on square or tri- 
angular lattices. In the upper half plane of ^, there is 
no quantum Monte Carlo study on this model due to the 
minus-sign problem. Other calculations based on mean 
field theory and exact diagonalization showed that the 
phase 7r/4 < < 7r/2 is antiferro-quadrupolar ordered 
on the triangular or square lattice^^. 



II. METHODS 

The tensor renormalization group method recently 
developed is an accurate numerical method for study- 
ing the ground state of quantum lattice models in two 
dimensions^^^^. It does not have the minus-sign prob- 
lem encountered in the quantum Monte Carlo simulation 
and can be used to study the phase diagram in the whole 
parameter space. We assume that the ground state is 




Figure 1: (color online) The ground state phase diagram 
of the spin-1 bilinear-biquadratic Heisenberg model on the 
honeycomb lattice. The inner circle is the phase diagram 
in the classical limit, while the outer circle is for the cor- 
responding quantum spin model. FM, AFM, FQ, AFQ 
and PVBS stand for ferromagnetic, antiferromagnetic, ferro- 
quadrupolar, antiferro-quadrupolar, and plaquette valence 
bond solid phases, respectively. Od ~ O.IOtt. 



described by the following tensor-product wave function 
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where rrii is the eigenvalue of spin operator 6'f. 
A\..y.^. [rrii] is the third-order tensors defined on the 6 
sublattices, as shown in Fig. [2l The trace is to sum 
over all spin configurations and all virtual bond variables. 
This wave function satisfies the area law of entanglement 
entropy. It is an accurate representation of the ground 
state wave function. Its accuracy is determined by the 
bond dimension D. It approaches the exact result in the 
limit D ^ oo. 

The ground state wave function, or the local tensors 
A*, is determined by applying the projection operator 
exp (y—rH) to an arbitrary initial state |^) iteratively un- 
til it is converged. Since this model only contains nearest 
neighbor interactions, exp {—rH) can be divided into a 
sequence of local two-site operators approximately by the 
Trotter- Suzuki decomposition for a sufficiently small r. 
We apply the first order Trotter- Suzuki decomposition 
here. In our calculation, we start the projection with 
a relatively large r = 0.2 and then reduce it gradually 
to 10~^ until the wave function is converged. In order 
to find the true ground state and not being trapped in 
a local minimum, we start the projection from variety 



Figure 2: (color online) Diagrammatic representation of the 
tensor- network wave function on the honeycomb lattice. Ten- 
sor A* defined on each lattice site contains three virtual bond 
indices and one physical index. 

of possible magnetically ordered states or valence bond 
solid states. We choose the converged state which has 
the lowest energy as the ground state wave function. A 
detailed introduction to this method can be found from 
Refs. [24.25]. This method is a fast and accurate way to 
get the ground state wave function. 

After obtaining the ground state wave function |^), we 
can evaluate the expectation value of physical variable O 



By contracting the physical indices, both (^| O |^) and 
(^1^) can be also expressed as tensor network. The con- 
traction of tensors is achieved by computing the dom- 
inant eigenvector of the corresponding one dimensional 
transfer matrix using the infinite time-evolving block dec- 
imation (iTEBD) method^^ beyond unitary evolution. 
The iTEBD is also an iterative projection method and 
the truncation error does not accumulate during the it- 
eration. The largest eigenvector of the transfer matrix is 
represented by a matrix product state with bond dimen- 
sion X, which determines the accuracy of the expectation 
values. 

In our calculations, we found that the ground state 
energy is converged when the bond dimension D > 12, 
while the expectation values of physical variables become 
stable when the parameter x ^ 30 (see Fig. [3]). Thus, we 
choose D = 12 and x = 30 throughout the calculations. 



Figure 3: (color online) The ground state energy density as a 
function of D with different values of x at ^ = 0. The curves 
of X = 30 and x = 40 are almost on top of each other. 
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Figure 4: (color online) The ground state energy density as a 
function of 0. 

and second derivative. The first derivative is calculated 
with Hellmann-Feynman theorem, which is more accu- 
rate than numerical differentiation from the ground state 
energy. We find that there are four phase transitions, 
located at ^ = — 37r/4, ±7r/2 and Od ~ 0.197r, respec- 
tively. Among them, 6 = — 37r/4 and ±7r/2 are first order 
transitions. The transition at Od is a second order one. 
This transition point is shifted below the classical value 
^ = 7r/4, which will be discussed later on. 

To clarify the phase diagram, we calculate various or- 
der parameters, i.e. the magnetization 



III. RESULTS 



The ground state energy shows that most part of the staggered magnetization 
the quantum phase diagram matches with the classi- 
cal phase diagram. Fig. [4] displays the 6>-dependence _ y^^_ji^^^^^^ (j^ 
of the ground state energy and Fig. [5] displays its first ■ 
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Figure 5: (color online) The first and second (the inset) 
derivative of the ground state energy with respect to 0. 
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Figure 6: (color online) Uniform (red line) and staggered (blue 
line) magnetization per site as a function of 6. 



the ferro-quadrupolar moment 

Q''=^{Sl)-2/3, (8) 

i 

and the antiferro-quadrupolar moment 

i 

in the four phases, respectively. Fig. [6] shows the - 
dependence of and M|. The ground state is found 
to have FM long range order for 7r/2 < 6 < 57r/4, and 
AF long range order for — 7r/2 < < Od- In these 
phases, the quadruple moment Q^^ is finite. In the re- 
gion — 37r/4 < 6 < — 7r/2, both magnetization and stag- 
gered magnetization vanish, however, the quadrupolar 
moment is finite, shown in Fig. [71 It corresponds to 
a ferro-quadrupolar phase, in agreement with both the 
semiclassical^^'^^ and quantum Monte Carlc^^ results. 
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Figure 7: (color online) Uniform (red line) and staggered (blue 
line) spin quadruple moment per site as a function of 0. 




(a) (b) (c) 

Figure 8: (color online) Pictorial representation of three pos- 
sible VBS patterns considered in calculation, (a) plaquette; 
(b) columnar; (c) staggered. The red thicker bonds represent 
stronger correlation while the black thinner bonds represent 
weaker correlation. 

In the staggered magnetization curve a sharp jump 
shows at Oc = —7t/2. This feature was also observed 
in the quantum Monte Carlo calculation on a square 
lattice^^. But the quadrupolar moment is finite and 
changes continuously at this point. So a first order 
phase transition occurs, in consistence with the conclu- 
sion drawn from the first derivative of the ground state 
energy. 

As expected, the staggered quadrupolar moment Q^^ 
vanishes in the FM, AF and ferro-quadrupolar phases. A 
surprising result is that this moment also vanishes in the 
classical staggered quadrupolar phase 7r/4 < 6 < 7r/2, 
i.e. the quantum fluctuation suppresses completely the 
staggered quadrupolar order, different from the previ- 
ous studies on the triangular or square lattices^^. More 
interestingly, the critical point has been shifted by the 
quantum fluctuation from 7r/4 to about O.IQtt, which ex- 
cludes the SU(3) AF Heisenberg spin-1 model from any 
long-range magnetic order. 

To further characterize the phase for 6d < < 7r/2, 
we have performed a thorough exploration of three pos- 
sible VBS patterns (Fig. [8]) on the honeycomb lattice. It 
has been checked that whatever VBS patterns we start 
with, it always converges into the plaquette order phase 
(Fig. [81(a)) under renormalization group flow. Hence the 
ground state energy of the plaquette VBS phase is the 
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Figure 9: (color online) The plaquette order parameter as a 
function of 0. 

lowest. This plaquette order phase explicitly breaks the 
lattice translation symmetry, but not the spin SU(2) sym- 
metry. A naive picture of this plaquette order is that, in 
order to minimize the ground state energy, the spins on 
one third of the minimal hexagons of the honeycomb lat- 
tice form the VBS phase, like the Haldane gapped phase 
in one dimension. 

In order to detect the plaquette order, we calculate the 
plaquette order parameter defined as 

E(i,j)Gred(^^ • ^i) ^ ^^^^ 

'^^{i,j)eblacki^i ' ^j) 

where (i, j) G red {black) means the two nearest neighbor 
spins connected by red (black) bond of Fig. [8fa). Fig. [9] 
shows the ^-dependence of the plaquette order parame- 
ters. Both the plaquette and AF orders vanish simulta- 
neously and continuously at the critical point Od. This 
observation suggests that this plaquette- AF transition is 
in fact a second order transition, in consistence with the 
conclusion drawn from the first and second derivative of 
the ground state energy with respect to in Fig. [5l But 
we are still unable to rule out the possibility of a very 
weak first order transition, partly due to the finite bond 
dimension D. 

However, around the transition point between anti- 
ferromagnetic and plaquette valence bond solid phases, 
we have evaluated the ground state wave function us- 
ing the cluster update approach proposed by Wang et 
al.^^. The cluster update considers long range entangle- 
ment by computing larger block size. The accuracy of the 
computation near a second order phase transition can be 
improved by using relatively small cluster size (6 sites). 



IV. DISCUSSION AND CONCLUSION 



A second order transition between AF and VBS order 
was originally predicted theoretically as the deconfined 
quantum critical point. This theory is based on the ob- 
servation that the topological defect (Skyrmion) of the 
AF order parameter carries a finite lattice momentum, 
thus after the AF order is suppressed by the Skyrmion 
proliferation, the system automatically enters the VBS 
order. The previous studies on deconfined criticality were 
focused on spin- 1/2 systems only, and this theory has 
gained strong numerical evidence from quantum Monte 
Carlo simulation on a spin- 1/2 model on the square lat- 
tice with both two-body and four-body interactions^ —. 
Our present result actually gives rise to a possible decon- 
fined quantum critical point in the spin-1 systems. Using 
the techniques in Ref. (33l35| . we can show that for spin-1 
systems on the honeycomb lattice the momenta carried 
by the Skyrmion will precisely lead to the plaquette order 
pattern after the Skyrmion proliferation^ 

Due to the critical point between AF and plaquette 
order at Od^ our numerical result also implies that the 
ground state of the SU(3) AF Heisenberg model on the 
honeycomb lattice has a plaquette order. This result con- 
curs with the recent studies on the SU(N) Heisenberg 
model^'', which suggested that the SU(N) spins tend to 
form block singlets that are commensurate with the lat- 
tice. 

To summarize, the ground state phase diagram of the 
quantum spin-1 bilinear-biquadratic Heisenberg model 
on a honeycomb lattice has been determined precisely. 
Besides the ferromagnetic, antiferromagnetic and ferro- 
quadrupolar phases, a plaquette order phase is found in 
the region of 6>(^ < 6> < 7r/2 for the first time, where the 
classical AF or staggered quadrupolar order is completely 
suppressed by quantum fluctuations. The quantum phase 
transition between AF and the plaquette order phase is 
found to be either a direct second order or a very weak 
first order transition. This is a possible candidate of a 
deconfined quantum critical point in a quantum spin-1 
system. Further investigation on the critical properties 
around this point is desired. 
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